Coherent storage of photoexcited triplet states using ^^Si nuclear spins in silicon 
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Pulsed electron paramagnetic resonance spectroscopy of the photoexcited, metastable triplet state 
of the oxygen- vacancy center in silicon reveals that the lifetime of the ms=dil sub- levels differ 
significantly from that of the rris—O state. We exploit this significant difference in decay rates to 
the ground singlet state to achieve nearly ^100% electron spin polarization within the triplet. We 
further demonstrate the transfer of a coherent state of the triplet electron spin to, and from, a 
hyperfine-coupled, nearest-neighbor ^^Si nuclear spin. We measure the coherence time of the ^^Si 
nuclear spin employed in this operation and find it to be unaffected by the presence of the triplet 
electron spin and equal to the bulk value measured by nuclear magnetic resonance. 

PACS numbers: 76.30.-v, 71.55.Cn, 76.70.Dx, 03.67.Lx 



Nuclear spins in solids are promising candidates for 
quantum bits (qubits) as their weak coupling to the envi- 
ronment often leads to very long spin coherence times [l]- 
[4|. However, performing fast manipulation and control- 
ling interaction between nuclear spin qubits is often more 
challenging than in other, more engineered, quantum sys- 
tems [5]-[7]. The use of an optically driven mediator spin 
has been suggested as a way to control coupling between 
donor electron spins in silicon: the donor spins exhibit 
weak direct coupling, but mutually couple through the 
optically excited state of the mediator [8]. Such ideas 
could similarly be applied to couple nuclear spins, and, if 
the mediator spin is a photo-excited triplet with a spin- 
zero single ground state, it would have the added advan- 
tage that it avoids long-term impact on the nuclear spin 
coherence [9UTT]. 

Photoexcited triplets are optically-generated electron 
spins {S = 1) which often exhibit large (positive or neg- 
ative) spin polarization, thanks to preferential popula- 
tion of each of the triplet sub-levels following intersystem 
crossing and/or the differing decay rates of these sub- 
levels to the ground singlet state [I?, "iT. Nuclear spins, 
in contrast, have weak thermal spin polarization at ex- 
perimentally accessible conditions, due to its small mag- 
netic moment. Highly polarized electron spin triplets can 
be used to polarize surrounding nuclear spins, through 
continuous wave microwave illumination (under processes 
termed dynamic nuclear polarization) [Ml [15] , or using 
microwave pulses [16 . Triplet states can also be used to 
mediate entanglement between mutually-coupled nuclear 
spins [3, on timescales much faster than their intrinsic 
dipolar coupling [17 . 

Oxygen- vacancy (0-V) complexes can be formed in sil- 
icon by electron beam or 7-ray irradiation of oxygen-rich 
silicon crystals |19] , and can be excited to the triplet 



state (termed an SLl center,) using illumination of above 
band gap light [20 . Magnetic resonance studies includ- 
ing electron paramagnetic resonance (EPR), electrically 
or optically detected magnetic resonance, spin depen- 
dent recombination and electrically-detected cross relax- 
ation [2QH25] have revealed that the SLl has orthorhom- 
bic symmetry with non-equilibrium triplet electron spin 
polarization and strong hyperfine coupling with the near- 
est neighbor ^^Si nuclear spins. The triplet spin polari- 
sation of SLl centers can be incoherently transferred to 
bulk ^^Si nuclear spins in the lattice by all-optical meth- 
ods [26 or dynamic nuclear polarization [27l [28] . Previ- 
ous electron spin echo studies on SLl performed at zero 
magnetic field have revealed that the populating rates 
of the triplet sub-levels are equal, and the spin polariza- 
tion arises instead from a difference in the decay rates to 
the singlet ground state [29]. In this Letter we use the 
high spin polarization of the triplet system and its strong 
coupling with the nearest neighbor ^^Si nuclear spins to 
demonstrate coherent state transfer between the electron 
and nuclear spin degrees of freedom, and examine the nu- 
clear spin coherence in the presence of the triplet. 

Cz-grown, single-crystal natural silicon (4.7% ^^Si, 
7 = 1/2) was exposed to 1 MeV e-beam irradiation (dose 
^ 10^^ cm~^) at room-temperature to form 0-V com- 
plexes (an interstitial oxygen already present in the sili- 
con traps a mono- vacancy generated due to the e-beam 
irradiation). Pulsed EPR measurements were carried out 
at X-band (9.72 GHz) on a Bruker Elexsys580 spectrom- 
eter equipped with a helium- flow cryostat. Photoexci- 
tation of the SLl was achieved using a 1064 nm pulsed 
NdiYAG laser (pulse width ^ 7 ns, 1 mJ/pulse) with a 
10 Hz repetition rate. 

Figure [1^ illustrates the SLl center in silicon under two 
representative orientations of the static magnetic field 
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FIG. 1: (color online), (a) Structures of the oxygen- vacancy 
centre in silicon, illustrating the SL1° and SL1^° orientations 
with respect to the externally applied magnetic field ^o- (b) 
The SLl triplet {S = 1) state is Zeeman split by Bq into levels 
T+, To, and T_. These states decay with different rates (A:+, 
/co, and k-, respectively) to the ground singlet {S = 0) state. 
The hyperfine coupling to the ^^Si (/ = 1/2) nuclear spin at 
lattice site i or j further splits the triplet states, (c) Electron 
spin echo-detected EPR spectrum of the SLl center at 12 K 
with Boll (110). The satellite peaks (red arrows) arise from 
hyperfine coupling to ^^Si. 



Bq. Under one orientation, termed SLl^, the magnetic 
field lies in the plane, marked (110), comprising the oxy- 
gen atom and two vacancy-trapping silicon atoms {i and 
j lattice sites). An alternative orientation (SLl^^) has 
the magnetic field in an orthogonal plane (110) with re- 
spect to the same center. Both planes are equivalent by 
symmetry for the crystal as a whole, so both orientations 
are visible in the EPR spectrum. 

Figure [TJd shows the triplet energy sub-levels in the 
presence of a static magnetic field (T+, Tq, and T_), each 
of which decays with a characteristic rate to the ground 
singlet state [24', ^30' . The hyperfine coupling with one of 
the two nearest-neighbor ^^Si nuclear spins (occupying i 
or j lattice site) further split the T±i sub- levels, while 
the To state {rris = 0) has no first-order hyperfine inter- 
action, thus the nuclear spin splitting in this sub-level is 
close to the Zeeman energy of 2^ Si. The EPR spectrum 
obtained by monitoring the electron spin echo intensity 



as a function of magnetic field at 12 K with Bq\\ (110) 
is shown in Fig. [l]^, labeled with electronic transitions 
identified by previous cw-EPR studies [20] . The satellite 
peaks accompanying each main peak are due to the hy- 
perfine interaction with the ^^Si nuclear spins situated at 
i or j lattice sites. The phase difference of the spin echo 
(i.e. dips or peaks) is indicative of the non-equilibrium 
polarization within the electron spin triplet. 

To investigate the origin of this non-equilibrium po- 
larization, we studied the decay kinetics of the triplet 
by measuring the electron spin echo at a variable time 
T after the optical excitation {hiy-T-7r/2-r-7T-r-echo)^ as 
shown in Fig. |2^. The zero echo intensity at T = in- 
dicates equal initial filling of the three triplet sub-levels 
upon creation of the triplet. The echo intensity proceeds 
to grow as T is increased. This can be attributed to a 
difference in the decay rates of the triplet sub-levels to 
the ground state, creating a build-up in spin polarisation 
(positive or negative) across the EPR transition being 
measured. 

Based on the simple decay model shown in Fig.[l]3, the 
population difference between a given pair of sub-levels 
follows a biexponential behavior where the two time con- 
stants represent the lifetimes of the two sub-levels in- 
volved in the EPR transition. The time constants ob- 
tained from biexponential fitting to the FD curves are 
given m Table [I]— the assignment of rates to particular 
energy levels is enabled by electron nuclear double reso- 
nance experiments described further below. As the tran- 
sition from the triplet to the ground singlet is determined 
by the amount of singlet admixture to the triplet via spin- 
orbit coupling, the lifetimes are expected to depend on 
the defect orientation with respect to the magnetic field. 
The composition of tripet levels T+^o,- can be expressed 
in terms of the zero-field eigenstates (T^^^^^^), and simi- 
larly the observed decay rates from these levels can be 
traced back to a corresponding mixture of zero-field de- 
cay rates (kx^y^z)^ as shown in Table 1. These values are 
in good agreement with times measured using zero-field 
EPR [29 . 

Our model assumes that there is negligible spin-lattice 
relaxation within the triplet sub- levels, and this is con- 
sistent with lack of temperature dependence we observe 
in the relaxation dynamics below 20 K. Neverthless, in 
order to probe the dynamics in more detail, we can in- 
troduce an additional inversion tt pulse to the sequence: 
{hu-Tx-TT-TY-TT /2-T-7r-T-echo). Figure [2]3 shows the 2D 
plot of the echo intensity for this sequence, as both Tx 
(the delay after the laser pulse) and Ty (the delay after 
the inversion pulse) are varied. The simulation of this 
experiment, based on the model described above, is in 
good agreement with the observed behavior, supporting 
our assumption that spin-lattice relaxation can be ne- 
glected. 

Based on the observed decay rates, we can extract 
both the polarization buildup (Fig. ^) and the triplet 
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FIG. 2: (color online), (a) Decay traces obtained by the flash 
delay {hv-T-n /2-T-Ti-T-echo) pulse sequence, used to extract 
the triplet decay rates shown in Table [l] (b) Experimental and 
(c) simulated 2D plots of the decay characteristics observed 
with the pulse sequence hiy-Tx-7v-TY-7r/2-T-n-T-echo. (d) The 
polarization build-up and (e) triplet population as a function 
of waiting time T after the laser pulse. 



population (Fig. |2p) for the two SLl orientations (SLl^ 
and SLl^^) as a function of time T after the laser pulse. 
The maximum electron polarization reaches > 99% after 
about 1.5 ms following the laser pulse. Below we investi- 
gate the coherent transfer of such well-prepared electron 
spin states to a neighboring ^^Si nuclear spin. 



TABLE I: Lifetime of triplet sub-levels, for two orientations 
of SLl centers, obtained from fitting to the flash delay curve. 



EPR transition 


Lifetime 


^ +1, SLl^ 
^ -1, SL1° 


(l/ko)^= 2000(4) /is, (l/k+i)^= 280(2) /is 
(l/ko)°= 1970(4) /is, (l/k_i)°= 330(2) fis 


-1 ^ 0, SLl^^ 
+ 1 ^ 0, SLl^o 


(l/k_i)^^= 960(2) /is, (l/ko)^^= 200(1) /is 
(l/k+i)90= 987(2) /is, (l/ko)^°= 205(1) /is 



= L6(3) ms~\ ky = 4.93(6) ms~\ ^ 0.50(4) ms~' 



For the following electron nuclear double resonance 
(ENDOR) experiments, we focus on four selected lev- 
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FIG. 3: (color online), (a) Davies-ENDOR spectrum illustrat- 
ing the hyperfine coupling between the triplet and nearest- 
neighbour ^^Si for the SL1^° center, (b) Rabi oscillation of 
the ^^Si nuclear spin, driven between the states |1) and |2), 
detected by monitoring the electron spin echo intensity. 



els (labeled |1) , |2) , |3) and |4) in Fig. [T]d) of SLl^^^ We 
studied the hyperfine coupling strength between the ^^Si 
nuclear spin and triplet electron spin using the Davies 
ENDOR pulse sequence (Fig.|3^). For SLl^^, Tq decays 
more quickly than T± (see Table thus the states |3) 
and 1 4) are mostly unpopulated in ~700 jas. A selective 
microwave 7r-pulse between |1) and |3) creates a polar- 
ization across the nuclear spin transitions, which can be 
driven using a radiofrequency (i^^f) pulse. The ENDOR 
signal is obtained by monitoring the electron spin echo 
on the |1):|3) transition as a function of Urf- Fig.|3^ shows 
the |1):|2) transition frequency dominated by the strong 
hyperfine interaction. Thus, with resonant rf pulses we 
can selectively address ^^Si nuclear spins at specific lat- 
tice sites i and j. These nuclear spins can be coherently 
manipulated with high fidelity, as illustrated by the Rabi 
oscillations in Fig. [3]3, in addition to being prepared and 
measured using the triplet electron spin. 

Using a sequence based on Davies ENDOR, it is pos- 
sible to coherently transfer a state of the electron spin to 
a coupled nuclear spin (see Ref [4 for full details). We 
apply this to transfer the highly polarised triplet elec- 
tron spin coherence to, and from, the nearest neighbour 
^^Si nuclear spin (see Figure [4^). The decay of the recov- 
ered spin coherence as a function of the storage time in 
the nuclear spin (2rn) is shown in Fig.|4]3, with an expo- 
nential decay of time-constant 0.9(1) ms. The measured 
decay is dominated by the relaxation of the T+ sub-level 
back to the ground singlet state (l//c+ = 987 /is), rather 
than ^^Si nuclear decoherence. This confirms our assign- 
ment of the decay rates shown in Table [ij By subtracting 
/c+ from the decay in Fig.|4]3 we can estimate T2 of ^^Si to 
be several milliseconds, however, it is possible to make a 
more accurate measurement as described below. 

Using the pulse sequence shown in Fig.[4]3, we shift the 
RF pulses within a fixed experimental window of 1 ms 
in order to remove the effect of the triplet relaxation and 
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FIG. 4: (color online), (a) The electron-spin coherence (EC) 
of the triplet is transferred to, and from, nuclear co- 

herence (NC) . (b) The retrieved electron spin echo intensity 
(blue) decays as a function of the ^^Si nuclear storage time, 
2Tn . A fit to an exponential decay (red) yields a time constant 
of 0.9(1) ms, apparently limited by the lifetime (0.99 ms) of 
the T+ state, (c) By fixing the experimental time window to 
^ 1 ms, the intrinsic coherence time of the nuclear spin can be 
measured, rf pulses were phase-cycled to confirm here that 
the measured electron spin echos arose solely from nuclear 
spin coherences. 



directly measure the ^^Si coherence time [31 . The se- 
quence is based on Davies ENDOR described above, but 
with the single rf tt pulse replaced with a nuclear Hahn 
echo sequence, whose delay time is swept. In the ab- 
sence of nuclear spin decoherence, the applied rf pulses 
form a net 27r rotation. In contrast, when the nuclear spin 
is fully decohered, the nuclear spin polarisation across the 
|1):|2) transition falls to zero. Fitting the data to an ex- 
ponential decay gives the nuclear coherence time of 5(1) 
ms. The bulk value for of ^^Si in natural silicon has 
been measured by NMR and found to be 5.6 ms, limited 
by ^^Si dipolar coupling [32 . Remarkably, the nuclear 
spin coherence appears unaffected by strong coupling to 
the triplet electron spin. 

In conclusion, we utilized the coupling between nuclear 
spin and photoexcited electron spin triplet in silicon to 
demonstrate the coherent storage and retrieval of triplet 
electron spin coherence in the ^^Si nuclear spin. This 
motivates further studies in how the nuclear spin state 
survives the decay to the ground singlet state, as well 



as the application of NMR pulse sequences to remove 
the effect of nuclear spin dipolar couplings to achieve 
coherence times of up to 25 seconds [3 . Furthermore, 
given the well-developed silicon isotope engineering ^31 
|34] it will be possible to investigate more than one nuclear 
spin strongly coupled to the single electron spin in the 
photoexcited triplet state and explore optical control of 
the interaction between the nuclear spins. 
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